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We report high-precision measurements of the Nusselt number Nu as a function of the 
Rayleigh number Ra in water-filled rectangular Rayleigh-Benard convection cells. The 
horizontal length L and width W of the cells are 50.0 cm and 15.0 cm, respectively, and 
the heights H = 49.9, 25.0, 12.5, 6.9, 3.5, and 2.4 cm, corresponding to the aspect ratios 
(T x ee L/H, T y ee W/H) = (1, 0.3), (2, 0.6), (4, 1.2), (7.3, 2.2), (14.3, 4.3), and (20.8, 6.3). 
The measurements were carried out over the Rayleigh number range 6 x 10 5 < Ra < 10 11 
and the Prandtl number range 5.2 < Pr < 7. Our results show that for rectangular 
geometry turbulent heat transport is independent of the cells' aspect ratios and hence 
is insensitive to the nature and structures of the large-scale mean flows of the system. 
Q-f This is slightly different from the observations in cylindrical cells where Nu is found to 

be in general a decreasing function of T, at least for T = 1 and larger. Such a difference 
is probably a manifestation of the finite plate conductivity effect. Corrections for the 
influence of the finite conductivity of the top and bottom plates are made to obtain the 
estimates of Nuoo for plates with perfect conductivity. The local scaling exponents /3; of 
Nuoo ~ Ra@ l are calculated and found to increase from 0.243 at Ra ~ 9 x 10 5 to 0.327 
at Ra ~ 4 x 10 10 . 



1. Introduction 

Convection is ubiquitous in nature and in our everyday life. It can be found in the 
stars and planets, in the Earth's mantle and outer core, in the oceans and atmosphere, 
as well as in heat transport and mass mixing in many engineering applications. The 
paradigmatic example for natural convection is the Rayleigh-Benard (RB) convection of 
an enclosed fluid layer between the colder top and the warmer bottom plates (Ahlers, 
Grossmann & Lohse 20096; Lohse & Xia 2010). A key issue in the study of turbulent RB 
convection is to understand how heat is transported upwards across the fluid layer by 
convective flows. The global heat transport by convection is usually expressed in terms 
of the Nusselt number, namely, 

(i-i) 

where Q is the heat current density across a fluid layer of thermal conductivity A/ 
with height H and with an applied temperature difference A. The dynamics of the 
system is determined by the geometrical configuration of the convection cell and by two 
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dimensionless control parameters, i.e. the Rayleigh number and the Prandtl number, 

aqAH 3 , „ v ,„ , 

Ra=— and Pr = -. 1.2) 

VK K 

Here, g is the acceleration due to gravity and a, v, and k are the isobaric thermal expan- 
sion coefficient, the kinematic viscosity, and the thermal diffusivity of the working fluid, 
respectively. The cell geometry is usually described in terms of one or more aspect ratios, 
such as r = D/H for a cylindrical cell of inner diameter D and (T x = L/H,T y = W/H) 
for a rectangular cell of horizontal length L and width W. The Ra- and Pr-dependence 
of Nu for various working fluids and cell geometries have been studied, both experimen- 
tally and numerically, in great detail for many years (Castaing et al. 1989; Kerr 1996 
Chavanne et al. 1997, 2001; Du & Tong 2000; Kerr & Herring 2000; Niemela et al. 2000 
Ahlers & Xu 2001; Roche et al. 2002, 2005; Xia, Lam & Zhou 2002; Verzicco & Camussi 
2003; Niemela & Sreenivasan 2003; Shishkina & Wagner 2007; Ahlers, Funfschilling k Bo- 
denschatz 2009a; Funfschilling, Bodenschatz & Ahlers 2009; Song & Tong 2010; Stevens, 
Verzicco & Lohse 2010; Stevens, Lohse & Verzicco 2011a; Silano, Sreenivasan & Verzicco 
2010; He et al. 2012). In addition, various theor etical models have been advanced to pre- 
dict the behavior of convective heat transport (ICastaing et al. 1989t Shraiman fc Siggia 



1990t iGrossmann fc Lohsel[2000l l200ll 120031 |2004 l201ltbubrullell200lU2002l ). For more 



detaile d elucidation of the problem, we refer interested readers to the recent review pa- 
per bv I Ahlers eial\ (|20096l ). O n the other hand, there are fewer measurements focusing 



on the aspect-ratio-dependence (Xu, Bajaj & Ahlers 2000; Fleischer & Goldstein 2002; 
Cheung 2004; Nikolaenko et al. 2005; Funfschilling et al. 2005; Sun et al. 2005a; Niemela 
& Sreenivasan 2006; Roche et al. 2010), and those measurements were all performed in 
containers with a cylindrical geometry. The objective of the present experimental investi- 
gation is to fill this gap by making high-precision measurements of Nu over a wide range 
of the aspect ratio in convection cells with a rectangular geometry. 

A lateral sidewall is indispensable for any convection experiment in the laboratory. The 
interaction between the sidewall and fluids would change the velocity and temperature 
distributions in the cell, and in turn change the flow structures of the system. Indeed, 
previous experimental studies for both cylindrical (du Puits, Resagk & Thess 2007) and 
rectangular (Xia, Sun & Cheung 2008) cells have shown that with increasing the aspect 
ratio the large-scale circulation (LSC) departs from a single-roll structure and becomes a 
multi-roll pattern. Thus, the T-dependence of Nu may reflect the influence of flow struc- 
tures on heat transport characteristics via the influence on the boundary layers. Mea- 
surements in cylindrical samples using water as working fluid (Pr » 4.3) revealed that 



for r < 6 Nu decreases , albeit only a few percent, with increasing T (jFunfschilling et al 



20051 : Sun et all 12005 ah . This suggests that the global heat transport properties of the 



RB system are not very sensitive to the flow structures for such parameter ranges. How- 
ever, the situation is very different for lower Pr and for the two-dimensional (2D) case. 
Using 3D direct numerical simulation (DNS), Bailon-Cuba, Emran & Schumacher (2010) 
found for Pr = 0.7 that the minimum of Nu occurs at the aspect ratio where the LSC 
undergoes a transition from a single-roll to a double-roll structure and the variations in 
Nu for different r are signi fi cant and can yield 11.3% for Ra = 10 8 . For 2D steady-state 
calculations, Idling & Tffl] (|2006h obtained a power-law of Nu ~ T 1 for T ^ 3. In 2D 
numerical RB flow, van der Poel, Stevens & Lohse (2011) identified different turbulent 
states for both Pr = 0.7 and 4.3, corresponding to different roll structures and associated 
with different overall heat transfers. Transitions among these states thus lead to jumps 
and sharp transitions in Nu(T). B y connecting the st r uctur es of iVti(r) to the way the 



flow organizes itself in the sample, Ivan der Poel et all (|201ll ) explained why the aspect- 
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Figure 1. Schematic diagram of front view of a rectangular cell. A: the Plexiglas cover; B: the 
copper top plate; Cj and C»' (i = 1,2,3,4): nozzles connecting the channels to the refrigerated 
circulators; D: thermistors; E: the Plexiglas sidewall; F: the copper bottom plate; G: the copper 
cover for the bottom plate, Hi: nozzle for transferring fluid into the cell; H2: nozzle for letting 
air out of the cell. 



ratio dependence of Nu is more pronounced for small Pr. Compared with the 3D cases, 
we note that the 2D simulations show larger variations in Nu(T) for both Pr = 0.7 and 
4.3, which may be explained by the different flow structures formed in the 2D and 3D 
cases. 

Different geometrical shapes represent different symmetries, and hence may lead to dif- 
ferent features of the flow and heat transfer. This prompts us to study the T-dependence 
of Nu in a non-cylindrical system. In the present study, we choose rectangle as the shape 
of the cells, which has been widely used in the past (Xia, Sun & Zhou 2003; Gasteuil et 
al. 2007; Maystrenko, Resagk & Thess 2007; Zhou k Xia 2008, 2010a,6; Zhou et al. 2007, 
2010). It was found that the convective flows in cells with such geometry share some simi 



lar dy namics as those in 2D RB system, such as the reversals of the LSC (Su givama et > 
l201ffl . 

The remainder of this paper is organized as follows. We give detailed descriptions of 
the experimental apparatus and conditions in §2. Experimental results are presented and 
analyzed in §3, which is divided into two parts. In §3.1, we compare the measured Nu 
for different aspect ratios and with those obtained in cylindrical cells. In §3.2 we consider 
the finite conductivity corrections of the top and bottom plates and estimate Nuoa for 
perfectly conducting plates. The scaling behaviors of Nu^ are & l so discussed in §3.2. We 
summarize our findings and conclude in §4. 



2. Experimental apparatus and methods 

Figure Q] is a schematic drawing of the front view of our apparatus and the drawing 
corresponds to (T x ,T y ) = (1,0.3). The sidewall of the cell, indicated as E in the figure, 
is composed by four transparent Plexiglas plates of 1.2 cm in thickness. The inner length 
L and inner width W of the cell are 50 cm and 15 cm, respectively. Six sidewalk of 
heights H = 49.9, 25.0, 12.5, 6.9, 3.5, and 2.4 cm were used in the experiment. The 
corresponding aspect ratios are (T x = L/H,T y = W/H) = (1,0.3), (2,0.6), (4,1.2), 
(7.3,2.2), (14.3,4.3), and (20.8,6.3), respectively. As T y is proportional to T x , in the 
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Figure 2. The normalized horizontal temperature differences, (T, — T m )/A, for the top (a, c) 
and bottom (6, d) plates and for T x = 1 (a, b) and 7.3 (c, d). The index i (= 1, 2, 3, 4, 5) of the 
thermistor is listed in figure [T] (see the top plate B and the bottom plate F). 



remainder of the paper we use only T x to indicate the cell's aspect ratio for ease of 
presentation. The top (B) and bottom (F) plates are made of pure copper of 56 cm 
in length and 21 cm in width and their fluid-contact surfaces are electroplated with a 
thin layer of nickel to prevent the oxidation by water. The thickness of the top plate is 
3.5 cm and that of the bottom plate is 1.5 cm. Silicon O-rings are placed between the 
copper plates and the sidewall plates to avoid fluid leakage. Eight stainless steel posts 
(not shown) hold the top and bottom plates together. They are insulated from the plates 
by Teflon sleeves and washers. Four parallel channels (not shown) of 1.5 cm in width 
and 2 cm in depth are machined into the top plate and the separation between adjacent 
channels is 1.5 cm. The channels start and end, respectively, at the two diagonal ends 
of the long edges. A silicon rubber sheet (not shown) and a Plexiglas plate (A) of 1.4 
cm in thickness are fixed on the top to form the cover and also to prevent interflow 
between the adjacent channels. At two ends of the i th channel (i — 1,2,3,4), there are 
two nozzles (Ci and CV) located, through which the channel is connected to a separate 
refrigerated circulator (Polyscience 9712) that has a temperature stability of 0.01 °C. 
The channels and the circulators are connected such that the incoming cooler fluid and 
the outgoing warmer fluid in adjacent channels always flow in opposite directions. To 
provide constant and uniform heating, two rectangular Kapton film heaters of 25 cm 
in length and 15 cm in width are sandwiched between two copper plates (F and G) 
and are connected in parallel to a d.c. power supply (SGI 330X15D) with 99.99% long- 
term stability. Therefore, the experiments were conducted under constant heating of the 
bottom plate while maintaining a constant tem perature at the top plate. Not e that recent 
high- resolution 2D ([Johnston fc Doering||2009l ) and 3D ( Stevens et aLll2011a ) simulations 
have revealed that turbulent thermal convection with boundary conditions of constant 
temperature and constant heat flux display identical heat transport at sufficient high 
Rayleigh numbers. 

Degassed water was used as the convecting fluid and the cell was leveled to better 
than 0.1°. During the measurements the entire cell was wrapped with several layers of 
Styrofoam. The temperature of each conducting plate was measured by five thermistors 
(D), which are embedded uniformly beneath the fluid-contact surface of the respective 
plate. When calculating the temperature difference A between the bottom and top plates, 
a correction was made for the temperature change between the thermistor position and 
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the fluid-contact surface. In each measurement after Ra was changed it took about 4 to 8 
hours for the system to reach the steady state and we waited for at least 12 hours to start 
the measurements. A typical measurement lasted over 12 hours and more than 24 hours 
for low-A experiments (A < 4 °C). No long-term drift of the mean temperature in the 
plates was observed over the duration of the measurement and the standard deviations 
were less than 0.5% of A for all measurements. 

To see the temperature uniformity of each conducting plate, we plot in figure [2] the 
normalized temperature variation, (Tj— T m )/A, for both the top and bottom plates. Here, 
Ti (i = 1,2,3,4,5) is the time-averaged temperature measured by the i th thermistor in 
a given plate and T m is the mean value of all the five thermistors in the same plate. It 
is found that the = 1 cell has the largest variation of plate temperature and we plot 
the r x = 1 results in figures ^ a) and (6). The temperature variation for the other five 
T x are similar and we choose T x = 7.3 as an example and plot the results in figures [2jc) 
and (d). One sees that (T$ — T m )/A is at most 4% for = 1 and is less than 2% for 
all I^ = 7.3 measurements. We thus estimate that a systematic error of the order of 1% 
could be introduced into the measured Nu by this horizontal temperature inhomogeneity 
of the conducting plates. Several thermistors were placed around the cell sidewall and 
around the bottom plate to monitor the environment temperature, based on which heat 
leakages to the environment and conduction by the posts and Plexiglas sidewall were 
calculated. The relative leakages to the total heat current were found to become more 
significant with decreasing A and thus we kept heat leakages from all sources to be less 
than 7% of the total applied heat current by working with sufficiently large A. We found 
that the largest source of leakage, especially for the small aspect-ratio cells, is through 
the cell sidewall which is insulated by multi-layers of Styrofoam. The errors in calculating 
the leaks come mainly from uncertainties in the thermal conductivities of the material 
involved, which are estimated to be less than 20%. This translates into an uncertainty of 
less than 1.5% in the results of Nu. 



3. Results and discussion 

3.1. Nu vs. Ra for different T x 

The measured Nu with corresponding values of A, Ra, and Pr are given in tables [1] 
and [2] The T x = 1 measurements were made at mean temperatures of 20° C and 30° C, 
corresponding to Pr — 7.01 and 5.45, respectively, and the measurements for the other 
five values of were conducted at 31.3°C, corresponding to Pr = 5.25. Previous stud- 
ies have revealed that with increasing Prandtl number Nu first increases, reaches its 
maximum value at around Pr ss 4 (depends on Ra), and then decreases slightly or re- 



mains independent of Pr (lAhlers fc XvJl2001t iGrossmann fc Lohsdl200ll : Kia et «Z.ll2002 ; 
Silano et al. 2010t Stevens et all 2011 ah . Therefore, within the present Pr range, Nu is 
expected to depend very weakly on Pr. Indeed, as shown in figures [3] and [5l no signifi- 
cant differences are observed between the Pr = 5.45 (circles) and 7.01 (right-triangles) 
results. Another source of uncertainty in the measured Nu co uld be the non-Bou s sinesq 
effects, as some of our measurements were made at larger A. iFunfschilling et all (|2005h 
argued that the applied temperature difference A should be limited to < 15°C to strictly 
conform to the Boussinesq conditions. However, it can be seen from figure [5] and tables 
H] and H that some of our data have A much larger than 15°C. As we shall see below, 
the measured large-A data show the same trend as those of small A. This suggests that 
some of our data being n ot strictly Boussine sq will not change the main conclusions of 
the present work. Indeed, Ahlers et all (|2006h have shown that for water as the working 
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Figure 3. (color online). Compensated Nu/Ra 1 ^ 3 on a linear scale vs. Ra on a logarithmic scale. 
V, = 20.8; O, = 14.3; A, T x = 7.3; □, T x = 4; <, F x = 2;>,T X = 1 and Pr « 7.01; 
0, = 1 and Pr ~ 5.45. Note that the data are as measured, without the correction for the 
finite plate conductivity. 



Table 1. Experimental parameters and results for the F x = 1 cell. H ere, Nuoo is calcula ted 
using Eq. (|3.ip and the parameters a = 0.275 and 6 = 0.39 obtained by iBrown" et al\ (2005) m 
cylindrical cells of two sets of plates of different thermal conductivities (Cu and Al). It is also 
seen that aA < 0.01 for all cases, which is considered to be sufficiently small. Note that two 
points are listed per line and the data are listed in chronological order. 
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fluid the non-Boussinesq effects could only slightly reduce the measured Nu by at most 
1.4% for A = 40 K and thus Nu is rather insensitive against even significant deviations 
from the Boussinesq conditions. Hence, for completeness all data are listed in the tables 
and are plotted in the figures. 

Figure [3] shows Nu/Ra 1 / 3 as a function of Ra for all six T x . It is seen that data points 
for T x — 20.8 (down-triangles) collapse well on top of those for = 14.3 (diamonds) 
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Table 2. Experimental results for the T x = 2, 4, 7.3, 14.3, and 20.8 cells. Note that two points 



are listed per line and the data are listed in chronological order. 
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(within their overlap Ra— range), which in turn collapse well on top of those for T x = 7.3 
(up-triangles). This implies that all sets of data can be described by a single curve over 
such a wide range of r x , i.e., no significant IVdependence is observed. As discussed in §1, 
Independence of Nu essentially reflects the influence of flow structures on heat transport 
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characteristics. Indeed, using particle image velocimetry (PIV), IXia et all (|2008l ) have 
shown in rectangular cells that the number of the convection rolls depends systematically 
on the aspect ratio of the system: only one convection roll is observed in the r x = 1 and 
2 cells and the LSC breaks into (horizontally arranging) multi-roll structure for T x larger 
than or equal to 4. Therefore, our present results suggest that for rectangular geometry 
turbulent heat transport is very insensitive to the nature and structures of the large- 
scale mean flows of the system. We note that our present r esults are different from those 



obtained in 2D numerical case 
Pr = 0.7 and 4.3. For the 2D simulation 
enable larger heat transfer than those with n 



van der Poel et ali 120111 ) , which were made for both 
the stable states with n rolls are found to 



1 rolls for vertically arranging LSC rolls. 
One possible reason for this difference may be attributed to different alignments of the 
LSC rolls, i.e. the aspect ratios of the 2D simulation vary between 0.4 and 1.25 and the 
LSC rolls are stacked vertica lly, whereas for the present geometry we have horizontally 
stacked rolls (|Xia et alj|2008h . 

In figure |H we compare the present results with those obtained in cylindrical cells and 
at Pr « 4.4. We mainly consider two recent data se ts: one is f rom Funfschilling et al. 
(2005) (referred to as FBNA) and the other is from Sun et al. ( 2005al ) (referred to as 
SRSX). As our measured Nu is independent of the cells' aspect ratio, here we do not 
distinguish our data for different r^. Note that all data in figure 4 are as measured, 
without the correction for the finite plate conductivity. For the FBNA data, one sees 
that the r = 6 data of FBNA are in excellent agreement with ours, but their small-r 
data are a few percent larger. Nevertheless, we note that the r = 1 data of FBNA (down- 
triangles) show similar trends to ours for Ra < 10 9 and Ra > 10 10 and thus the small 
difference between the two sets of data may be attributed to different system errors. The 
SRSX data and ours are much closer. One sees that parts of the L = 5 and 10 data of 
SRSX agree well with ours and the others lie slightly below our measured Nu. 

One noticeable difference between the present results and those of FBNA and SRSX 
is worthy of note. For both FBNA and SRSX the larger-r results lie consistently below 
those of smaller ones, i.e. Nu is generally smaller for larger L. This is true even for the 
largest-r (i.e. T = 10 and 20) data of SRSX (see figure U]). Whereas, our measured Nu 
for all six values of T x fall into a single curve, i.e. Nu is essentially independent of T x 
for our results. To understand such a difference, we note that all data plotted in figure [4] 
have n ot been correct ed for the influence of the finite conductivity of the top and bottom 
plates (|Verziccdl2004l ). In a cylindrical cell of T = 0.5, Sun, Xi & Xia (20056) argued that 
because of the finite conductivity and finite heat capacity of the plates the azimuthal 
sweeping of the circulation plane of the LSC would make heat transfer more efficient than 
the case when the LSC is locked in a particular orientation. Moreover, Xi fc Xia ( 2008h 
studied the azimuthal motion of the LSC in cylindrical cells of T = 0.5, 1.0, and 2.3 and 
their results showed that the LSC's azimuthal motion is more confined in larger-r cells. 
If we generalize the above two findings to large T and taken them together, a possible 
scenario can be achieved: the larger-r cell confines the azimuthal sweeping motion of the 
LSC which in turn reduces the measured Nu. This scenario could be valid for cylindrical 
cells due to their azimuthal symmetry. But for rectangular cells, this is not expected to 
work because th e rectangular geometry has already locked the orientation of the LSC 
( Xia et aLll2003l ). Therefore, our present T^-independent results are consistent with the 
lack of the azimuthal sweeping motion of the LSC in rectangular cells. What we should 



stress is that the large-r cell contains multi- roll structures of the LSC (|du Puits et al 



2007tlXia et a/.ll2008HBailon-cuba et aZ.|[2010h and thus the azimuthal motion of the LSC 



and its influence on heat transfer should be more complicated for large T. However, as 
we shall see in figure El when corrections for the finite conductivity of the plates are 
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Figure 4. (color online) . Comparison among Nu/Ra- 1/S fr om the present work (0)> from 
iFunfschilling et al\ (|2005l ) (y, V = 1 ; x, V = 2; *, V = 3; □, T = 6; for clarity their V = 1.5 data 
are not shown) and from ISun eFaU (|2005al ) (A, r = 2; O, F = 5; <3, F = 10; O, r = 20). As 
our results show that for rectangular geometry Nu/Ra~ 1 /' i is independent of the cells' aspect 
ratio, here we do not distinguish our data for different T x . Note that all data are as measured, 
without the correction for the finite plate conductivity. 
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Figure 5. (color online), (o) Compensated Nuoo /Ra 1 ^ 3 on a linear scale vs. Ra on a logarithmic 



scale: symbols as figure [3] Here, Nuqq i s calc ulated using Eq. (|3.1[) and the parameters a = 0.275 
and b = 0.39 obtained by Bro wn et al. (2005) in cylindrical cells. For comparison, the uncorrected 
data Nu/Ra 1/3 of figure O are replotted in (6). 



made, F-dependence of Nuoo becomes weaker for both the FBNA and SRSX data, which 
suggests that the finite plate conductivity effect is indeed a major factor for the observed 
difference of the behaviors of heat transfer in rectangular and cylindrical cells. 

3.2. Nuoo vs. Ra for different T x 
iBrown et al. ( 20051) suggested an empirical correction factor f{X) = 1 — exp[— (aX) b ], 



namely 

Nu = f{X)N Uoo = {1 - exp[-(aX) b ]}N Uoo (Ra, Pr), (3.1) 



10 



Q. Zhou et al. 



to obtain estimates of the ideal Nusselt number Nu^ for plates with perfect conductivity 
from the measured Nu. Here, X = Rf/R p is the ratio of the effective thermal resistance of 
the working fluid, Rt = H / (XfNu), to the thermal resistance of the plates, R p = e/X p , X p 
(= 401 W/m K) is the conductivity of plates (Cu), X f (A/ = 0.614 W/m K for Pr = 5.25 
and Xf = 0.589 W/m K for Pr — 7) is the conductivity of water, and e (= 1.5 cm) is the 
mean thickness of the top (the part of th e top plate below the cooling channels is used 
here) and bottom plates (|Verziccol[200l . To apply the relation (EH) to our measured 
Nu, one needs to determine the values of a and b. T he best way to do thi s is to use plates 
of different thermal conductivities, as was done by iBrown et al. ( 2005 ). who used two 
sets of plates made of Cu and Al. However, the lack of aluminum-plate measurements 
in the present study prevents us to determine the values of a and b. Alt ernatively, to 



estima te Nuoo, we use the parameters a = 0.275 and b = 0.39 obtained bv lBrown et al 
(2005) in cylindrical samples of 50 cm in diameter. 

Figure^a) shows the calculated Nu^/ Ra 1 / 3 as a function of Ra. For comparison, the 
uncorrected data Nu/Ra 1 / 3 of figure[3]are replotted in figure \M.b). The corrected data 
in figure E^a) seem to display a small aspect ratio dependence, e.g., near Ra — 2 x 10 s 
the T x = 7.3 data lie slightly above the T^. = 4 data and near Ra — 2 x 10 9 the = 4 
data lie slightly above the = 2 data. We note that this dependence is consistent with 
the FBNA data (see figure where data points for T — 6 lie slightly above those for 
T = 3, which in turn lie slightly above those for V = 2). Nevertheless, the differences 
of Nuoa for the two adjacent-r^ sets of data are only a few percent, which is smaller 
than or comparable with the experimental errors. Hence, we may also conclude that no 
significant T^-dependence of heat transfer is observed. 

We want to e mphasize tha t here we adopted the parameters of a — 0.275 and b = 0.39 
determined by | Brown et al\ (|2005l ) for cylindrical samples with a diameter of 50 cm. 
Brown et all (|2005l ) have shown experimentally that a and b are independent of the 



aspect ratio, but vary with the diameter of the sample. In the present study, we chose 
the rectangle as the geometry of the convection cells, which can possibly have an influence 
on the finite plate correction, as we have discussed at the end of § 3.1. However, without 
the measurement performed with aluminum plates it is difficult to assess whether the 
geometry would (significantly) influence the values of a and b. By using different sets of 
a and b to perform the finite conductivity correction, we estimate that the uncertainties 
of a and b may yield a few percent uncertainty on the obtained Nuoo , which is the same 
order of the difference between Nu and Nuoo- We should also stress that the best way 
to estimate iVuoo i s to use plates with different thermal conductivities, as was done by 



Brown et al\ (|2005l ) for cylindrical samples. Therefore, new measurements with Al plates 



are essential for concretely settling the problem and this will be the objective of future 
studies. 

In figure [6] we directly compare our Nu^ and those of FBNA and SRSX. It is seen 
that our data display the similar Ra- dependence as the r = 1 data of FBNA, and an 
excellent agreement between our data and those of SRSX can be found for nearly all the 
overlap Ra range. As the three sets of data were taken independently from the samples 
with very different geometries, this agreement is just remarkable. 

Finally, we studied the power-law relation between Nuoo and Ra, namely, 

Nu oc =A(T x )Ra^\ (3.2) 

Table [3] displays the fitted results for the power-law relation (|3.2[) at fixed aspect ratios. 
One sees that as increases in general the prefactor A increases and the scaling exponent 
/3 decreases. Here, T^-dependences of A and f3 essentially reflect their i?a-dependence as 
our results do not reveal significant aspect-ratio dependence of Nuoo. We thus do not 
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for the finite plate conductivity. 



Table 3. Fitted parameters from equations (|3.2[) . 

1 2 4 7.3 14.3 20.8 

A(T X ) 0.074 0.059 0.108 0.173 0.212 0.242 
/3(T X ) 0.324 0.335 0.307 0.282 0.271 0.262 
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Figure 7. (color online). Local scaling exponent /3j of Nuoo(Ra), determined from a power-law 
fit over a sliding window of half a decade, as a fun ction of Ra for the present data (solid circles) 
and for the F = 1 data of Funfschilling et al. (2005) (open triangles). Dashed line marks Pi = 1/3 
for reference. Note that both sets of j3i were obtained from the corrected data. 



distinguish the data for different and take them together to study the scaling behaviors 
of Nu 00 (Ra). The local scaling exponent fa is obtained by a power-law fit, Nuoo ~ Ra 13 ' , 
to the data for Nuoo(Ra) within a sliding window that covers half a decade of Ra. Figure 
[7] shows the results for fa as a function of Ra. It is seen that fa roughly increases linearly 
with log Ra from fa = 0.243 at Ra w 9 x 10 5 to fa = 0.342 at Ra « 3.6 x 10 9 . At small 
Ra, our results differ very much from the DNS results for cylindrical samples of unit 
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aspect ratio by Wagner, Shishkina & Wagner (2012), who found that ft increases again 
as Ra decreases below 2 x 10 7 and grows to about 0.30 at Ra = 10 6 [see figure 2 of 
Wagner et al. (2012)]. The apparent differences in ft may be explained by the different 
Prandtl number in the two studies, i.e., Wagner et al. (2012) carried out their simulations 
at Pr = 0.786, while our measurements were made at Pr re 5.45. For higher Ra, j3i drops 
slightly and fluctuates around 0.32. The FBNA results for ft are also displayed together 
with ours in figure [7] for comparison. For Ra 10 s the two data sets both increase with 
Ra with ours being a little scatter. A source of uncertainty for our data could be the 
non-Boussinesq effect, i.e. the FBNA data were obtained in the strictly Boussinesq range, 
while some of our data are beyond the Boussinesq range (i.e. A > 15°C). However, as 
have discussed in §3.1, some of our data being not strictly Boussinesq will not change 
our main conclusions. What is worthy of note is that around Ra re 10 10 our measured ft 
has a value that is close to the value of 1/3. The exponent /3 re 1/3 was obtained before 
in cylindrical cells o f T = 1 by FBNA at Ra re 7 x 10 10 (see open triangles in figure 
EJ) and of T = 4 by iNiemela fc Sreenivasanl (|2006l) for Ra > 10 10 . Here, our results in 
rectangular cells seem to be qualitatively consistent with these findings 



4. Conclusion 

In conclusion, our high-precision measurements of Nu in rectangular cells with (T x , T y ) 
varying from (1, 0.3) to (20.8, 6.3) show that Nu is independent of the aspect ratio. This 
is slightly different from the observations by both FBNA and SRSX in cylindrical cells 
where Nu is found to be in general a decreasing function of T, at least for r ~ 1 and 
larger. Such a difference may be attributed to different azimuthal dynamics of the large- 
scale circulation (LSC) and is probably a manifestation of the finite plate conductivity 
effect. To make finite conductivity corrections, an empirical correction factor f(X) = 
1 — exy\—{aX) b ], t ogether with the parameters a = 0.275 and b = 0.39 obtained by 



Brown et al\ ((20051 ) in cylindrical samples, were adopted to estimate Nu^ for plates 
with perfect conductivity from the measured Nu. The obtaiend Nu^ were found to be 
consistent with the FBNA and SRSX data measured in cylindrical samples to only a 
few percent. The scaling behaviors between Nuoo and Ra were studied for all six aspect 
ratios. The local scaling exponents ft of Nu^ ~ Ra@ l were calculated and found to 
increase with increasing Ra. Around Ra re 10 10 our measured ft has a value that is close 
to the value of 1/3. 

This work was supported by the Natural Science Foundation of China (Nos. 11161160554, 
11002085, and 11032007) and Shanghai Program for Innovative Research Team in Uni- 
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